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Abstract 

The contribution of the spin-orbit interaction in Hartree-Fock calculation for 
closed shell nuclei is studied. We obtain explicit expressions for the finite range 
spin-orbit force. New terms with respect to the traditional spin-orbit expressions 
are found. The importance of the finite-range is analyzed. Results obtained with 
spin-orbit terms taken from realistic interactions are presented. The effect of the 
spin-orbit isospin dependent terms is evaluated. 
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1 Introduction 



The success of the shell model in predicting the nuclear magic numbers is related to the 
presence of a strong spin-orbit term in the nuclear average potential. Few years after the 
formulation of the nuclear shell-model evidences of spin-orbit terms in the nuclear 
interaction were identified by analyzing polarized proton scattering data off complex 
nuclei 0. 

The strong spin-orbit term of the nuclear average potential should be generated by an 
analogous term present in the nucleon-nucleon interaction. However, the connection be- 
tween realistic nucleon-nucleon interactions, i.e. those built to reproduce the two nucleon 
scattering data and the deuteron properties, and the effective interactions, those used in 
nuclear structure effective theories, is still unclear. 

In the present article we present a study on the relationship between the spin-orbit 
terms used in realistic nucleon-nucleon interaction and the analogous terms used in Har- 
tree-Fock (HF) calculations. Our work has been done within the non relativistic frame- 
work where the spin-orbit terms can be easily isolated. In effect, in our HF calculations, 
we have used the explicit spin-orbit terms of the Argonne-Urbana nucleon-nucleon realis- 
tic potentials @, ^ § describing, within the non relativistic framework, nucleon-nucleon 
elastic scattering data up to energies of about 300 MeV. 

The spin-orbit interactions commonly used in HF calculations are of zero-range type 
and, in general, they are parametrized following the expressions proposed by Skyrme . 
Also the Gogny interaction which has a finite range for all the other channels, uses a 
Skyrme-like expression for the spin-orbit term. 

We explicitly develop the expressions for a finite range spin-orbit interaction to be 
used in HF calculations. As expected, in addition to the direct term these expressions 
produce a contribution also in the Fock-Dirac exchange term of the HF equations. This 
term is not present when zero-range interactions are used. Also in the direct term there 
are some new contributions with respect to the expression obtained with the Skyrme 
interaction. 

In the next paragraphs we present the detailed expressions of the HF equations when 
finite range spin-orbit terms are considered, then we discuss the importance of the finite- 
range and the effect of using spin-orbit terms taken from realistic nucleon-nucleon inter- 
actions. Finally we draw our conclusions. 

2 The formalism 

In ref. we made an explicit presentation of the HF formalism with finite-range interac- 
tions. In the present article we extend the formalism in order to treat also the spin-orbit 
terms. For this reason we recall here only those parts of the formalism involving the 
spin-orbit terms. 
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The effective interaction used in our calculations has the form: 



V{rur2) = E Vp{n,r2) 0^(1,2) (1) 
p=i 

where to the first 6 components used in ref. we have added the spin-orbit terms 
defined as 0\l, 2) = L • S and 0^(1, 2) = L ■ S ti • ra with L = (ri - ra) x (pi - pa) and 
S = 1(0-1 + 0-2). 

To take advantage of the spherical symmetry of the problem, we describe the single 
particle wave functions by separating the angular and the radial parts 0fc(r) = k{Q)uk{r)/r 
where the subindex k indicates all the quantum numbers necessary to identify the state 
and fl the two angular coordinates 6 and (p. 

The specific expression of the single particle wave functions allows us to reduce the 
HF equations into a set of differential equations of the type: 



Uk{r) + Uk{r)uk{r) - Wk{r) = ekUk{r), (2) 



2mk ydr^ 
where we have defined: 

Uk{r) = Y.j dr'<(r') j dVt j dfi' k\ny{Q.')V{\v - r'\)k{n)i{n') Ui{r'), (3) 

i 

and 

Wk{r) = Y.j dr'<(r') j dVt j dfi' F (fi)^(^]')^(|r - r'|)I(fi)A;(fi') Ui{r)uk{r'). (4) 



While the interaction depends from the relative distance between two nucleons, the HF 
equations depend upon the distance of the particles from the origin of the reference 
system. The implementation of finite range interactions in the HF equations requires 
the separation of the coordinate variables in the interaction. For the central and tensor 
channels this separation is done by considering the interaction in coordinate space as 
Fourier transform of the interaction expressed in momentum space (see ref. for details). 
For the spin-orbit channels (p = 7, 8) we use a different strategy consisting in expanding 
in multipoles the interaction: 

V,{t,2) = 47r 5: iv^(ri,r2) y£^,(n)lLM(f2). (5) 

LM ^ 

with L = V2L + 1. lYiom. the previous equation, making use of the orthogonality of the 
spherical harmonics, we obtain a close expression for the coefficients of the expansion: 

V£(ri,r2) = y dcos0i2Vp(ri2)PL(cosM, (6) 

In the previous equation we have indicated with the Legendre polynomials and with 
6*12 the angle between ri and Y2- 
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The details of the calculations of the spin-orbit matrix elements are given in the 
Appendix. We obtain for the direct term in the HF equations the following result: 



V'2 



jk Uk + 1) -h {k + 1) - 



where 



-'fc 



1, p=7 
2 4, p=8 



K{ry)+UUry)\, (7) 

(8) 



and the potentials and U^g are given by: 
Wg(r,r') : 

and 



vi{ry)-\-vi{ry) 

O T 



Ulsir,r') = [Vo^(r,r')-^^Vf(r,r') 



p=7, 



p = 7, 



(9) 

(10) 



The function ^ci^) used in the previous equations has been defined as: 

p(r) , p=7 



where p'^(r) and p^{r) are the proton and neutron densities such as p'^(r) + p^{r) = p{r). 
The other function used in eq. (|l3) has been defined as: 



pLs{r) , P=7 



(12) 



where pls is the nucleon spin-density, 



Ji{ji + l)-k{k + l)-- 



uAr] 



(13) 



and Pisi^) P'isi''^) analogous functions for protons and neutrons respectively. 

For the exchange terms of eq. (|^) we obtain: 



[Wk{r)] 



where 



p=7,8 



^ki 



iL a=l,5 •' 



p = 7 



(14) 



(15) 



The new five functions e have been defined as: 



'kiL 



a = 1 5 , 



(16) 
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with Ch\l,K) given by: 



i{h + k + L)%^{L,K) 

k h K 
-1 1 



kik + 1) 



k h K 
-1 1 



-J k Ik K 
1 

k k K\l 
0-1 1 



I 1 K L 
)\l -1 

I K L\ 
1 -1 j ' 



a 



a 



U Ik K\ 




K 




I / 


J 








a = 4,5. 



In the previous expressions we have used the following definitions: 

3k ji L \ 



Ik \ jk 



K 1 L 



1 _i 

2 2 " 



4 I JA: 
/■ i ?• 



/fc K 

Ji Jk 9 



L 1 K 

Gh{L,K) = ^{l, + l^ + L + l)J2aL + L' + l)L'' 

L' 

In the eqs. (|l^ we have used five new potentials: 



y 

1 K L' 
1 -1 



rk^iL',K). 



u*{r')Vl{r,r')ukir')ui{r) \ ^, a 



' 1, 


a 


= 1, 








r' 
f ■> 


a 


= 2, 


r 

[ r" 


a 


= 3, 



W,1?(r,r') 



r'M*(r')Mfc(r') V£(r,r') ^Mi(r), a = 4, 

'"WfelO^KlO V£(r,r')]Mi(r), a = 5. 
Like in ref. M the numerical solution of eq. (0) has been obtained iteratively using the 



plane wave expansion method of refs. ||10[- The center of mass motion has been considered 
in its simplest approximation, consisting in inserting the nucleon reduced mass in the 
hamiltonian. The single particle wave functions used to start the iterative procedure have 
been generated by a Saxon-Woods potential without spin-orbit and Coulomb terms. 
Therefore the starting wave functions for spin-orbit partners are the same. 
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3 Results 



In the same spirit of the work of ref. |^ we are more interested in investigating the vahdity 
of the commonly used approximations rather than proposing a new effective interaction 
to be used in HF calculations. This study has been conducted by adding different kinds 
of spin-orbit terms to a basic interaction composed by the four central terms of the force. 
These terms are described as a sum of two gaussians: 

Vp{r)=j2A,,exp{-hr^), (17) 

■t=i 

with p = 1,2,3,4. The parameters of this part of the interaction, which we call Bla, 
are compared in tab. |l] with the parameterization Bl of Brink and Boeker ||rT|. The 



small differences are due to the fact that we have considered the Coulomb interaction and 
therefore we had to readjust the parameters of the force in order to reproduce the binding 
energy of "^He. We added to the Bla interaction a finite range spin-orbit term of gaussian 
form: 

V,{ru) = A,expi-bjrl^). (18) 

The finite range effects have been investigated by comparing the results obtained with the 
above interaction with those produced by adding to the Bla force a zero range spin-orbit 
term of the form: 

Vriru) = Ar6\r,-r2). (19) 

The straightforward insertion of this expression in our formalism gives a contribution 
exactly equal to zero. The reason of this result can be traced back to the fact that we 
have developed our expressions using L = (ri — r2) x (pi — P2). To get results different 
from zero for a zero-range spin-orbit interaction we set to zero the quantity Vf(r, r') in 
eqs. (1^) and (p!oD , and after inserting eq. ([T9| ) we obtained: 

l^L,Lsir,r') = ^nl,sir')S{r-r'). (20) 

The calculations done with this approach are labelled as z. 

In addition to these effective interactions we have also used spin-orbit terms taken 
from microscopic forces: the Urbana V14 [Q, Argonne V14 ^j, and Argonne V18 
potentials. Our study has been restricted to the investigation of the doubly magic nuclei 
12c, 1^0, 40Ca, 48Ca and ^ospb. 



The finite range interaction (|T8D has been used to study the role played by the various 
terms of the spin-orbit potential. In a first set of calculations, labelled as c, only the Uc 
term of eq. (|^) has been used. This is the only spin-orbit term present in shell-model 
calculations. In another set of calculations, denoted as d, we considered the full direct 
term, and, finally, the results identified with so have been obtained with all the spin- 
orbit terms. These calculations have been done by changing every time the parameters 
of the force (|r^) to reproduce the 6.3 MeV splitting between the protons Ip levels in 
^^O. The parameter 67 was fixed to the arbitrary value of 1.2 fm"^ and the fit of the 



6 



splitting was obtained by changing A7. The values of A7 obtained in this way are, - 
108.75, -107.50, -97.86 MeV for the c, d, and so calculations respectively. The three 
interactions do not differ very much as it is shown in the panel I of the figure. This result 
indicates that the largest contribution from the spin-orbit force is coming from the Uc 
factor of the direct term. Since the other terms are small we have explored the possibility 
of avoiding their explicit calculation by simulating their effects with a readjustment of 
the force parameters. This is the reason why each type of calculation has been done with 
a different parametrization of the force, each of them reproducing the same empirical 
quantity. 

In tab. we compare the binding energies obtained with our calculations with the 



experimental ones |T2[. In spite of the fact that we handle with a non-linear problem, 
the spin-orbit terms acts on the binding energies as expected. The main contribution 
to the binding energy is obtained by the sum of the single particle energies. In nuclei 
where all the spin-orbit partners are occupied the spin-orbit term lowers the energy of 
the 1 + 1/2 level and increases that of the / — 1/2 level, in such a way that the contribution 
to the nuclear binding energy is almost zero. In table |^ this is observed by looking at 
the values of the energies of ^^O and ^°Ca which are practically the same, independently 
from the spin-orbit force used. Clearly those nuclei where not all the spin-orbit partners 
are occupied are sensitive to the spin-orbit force, since the single particle energy of the 
last occupied level is lowered. The effect is seen in ^^C, ^^Ca and ^''^Pb where the binding 
energy increases, in absolute value, the stronger the spin-orbit force is. 

The quantity most sensitive to the spin-orbit interaction is the energy splitting be- 
tween spin-orbit partners levels. The energy splittings calculated for the various nuclei 
under investigation with the interactions proposed are compared in tabs. ^ and ^ with the 
Skyrme III [|T3[ results and with the empirical values [Q. The experimental spectrum is 
more compressed than the theoretical one. This fact is well known ||15[, and it is related to 
the intrinsic limitations of the HF theory in the description of an interacting many-body 
system. 

As expected, the splittings increase with increasing value of /. The splittings obtained 
with the zero range interaction z become larger than those obtained with finite range 
interaction as the mass number of the nucleus increases. The results obtained with zero- 
range Skyrme interaction do not present this effect. In the Skyrme interaction there 
are velocity dependent terms generating spin-orbit like contributions which add to those 
produced by the genuine spin-orbit term. These velocity dependent terms simulate the 
effects of the finite range. We observe that the value of the splittings obtained with 
the Skyrme III interaction are comparable with those obtained with our finite range 
interactions. 

/^From the comparison of the results of the c and d columns of tabs. |^ and § we infer 
information on the role of the terms U^g in eq. (J^). The inclusion of Uf^g increases the 
splitting for all the nuclei considered but the magnitude of this increase is rather different 
for the various nuclei. We should not consider in our analysis the nucleus ^^O since it has 
been used to fit the interaction. We notice that the addition of Uf^g produces quite small 
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differences in the splitting of Ca and Pb nuclei while they are remarkable in C and 
^«Ca. 

These results can be understood by considering that lA^g is related to the nuclear spin 
density, eq. (|13|). If we assume that the radial wave functions u{r) are the same for spin- 
orbit partners levels, the contribution of these two levels to the spin density is exactly 
zero. In real calculations these wave functions are slightly different, but the contribution 
to the spin density remains small. This explains the small increase of the splitting in 
^°Ca and the relatively large modifications produced in ^^C and ^^Ca. One should remark 
that only the unoccupied levels contribute to the spin density. For this reason the effect 
of U^g is relatively large with respect to that of Ufj in ^^C and ^^Ca where the number 
of single particle levels is relatively small. In a heavy nucleus like ^osp]^ there are many 
levels contributing in U^; and the effects of U^g produced by a single level is relatively 
small. 

The contribution of the exchange term can be seen by comparing the results of the 
d and so columns. The variations with respect to the calculations done with only the 
direct terms can be as big as 10-15%, but not all of them have the same sign. It seems 
that for all the p states the splitting is reduced when the exchange term is considered, 
but it is increased in the /, g and h states. The situation for the d states is even more 
complicated, since the splitting is reduced for the Id states in ^ospb it has increased 
for all the other d states. The contribution of the exchange term cannot be taken into 
account in calculations with the direct term only by modifying the force parameters. 

The values of the splittings produced by the Urbana (U), and Argonne Vig interactions 
are comparable with those of our interactions, while the Argonne V14 (A14) generates 
smaller values. The radial dependence of the spin-orbit terms of these interactions are 
shown in the panel II of the figure. It is remarkable that the results of U and A18 are 
similar in spite of the large difference in the depth. The depth value of A14 is intermediate 
between those of the previous two forces, but its splittings are smaller. Th U and A18 
forces have similar range, while that of A14 is smaller. These facts indicate that our 
calculations are more sensitive to the range of the interaction than to its minimum value. 
In effect we recall that, in our calculations, a zero-range interaction does not produce any 
splitting. 

The calculations done with the microscopic interactions include both spin-orbit and 
spin-orbit isospin terms. In order to study the importance of the isospin part of the spin- 
orbit interaction we have repeated each calculation leaving out this terms. The differences 
of the results obtained with the full interaction and those without the isospin part are 
very small. In order to avoid a long list of numbers we give in tab. |^, for each nucleus 
under investigation, the minimum, the maximum and the average difference, in absolute 
value, between the calculated splittings. It appears clear the relatively small importance 
of this term of the interaction. This fact can be understood considering that the major 
contribution to the spin-orbit interaction is coming from the direct Uf; term. For the 
spin-orbit isospin term of the interaction, the case p = 8, the Uf^ term contains a function 
which is given by the difference between the proton and neutron density distributions, eq. 
(pTl). In all the nuclei we have considered this difference is small and particularly small 
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in those nuclei having the same number of protons and neutron (^^O and "^^Ca). In effect 
the maximum differences are larger in ^*Ca and ^osp^-, tj^an in ^^O and ^°Ca. 

4 Summary and Conclusions 

In this article we have presented a formalism to treat finite range spin-orbit interactions 
in HF calculations. The finite range of the interaction generates additional terms with 
respect to the usual shell model expression. One of these is the contribution to the 
exchange Fock-Dirac term in the HF equation (H). Also in the direct (Hartree) term of 
this equation there is a new part: the Uls piece of eq. (|^. The major goal of our work 
was the investigation of the effects produced by these new components. This has been 
done by adding different type of spin-orbit terms to a fixed interaction active only in the 
four central channels. We have used a spin-orbit interaction of a gaussian form whose 
parameters have been fixed to reproduce the energy splitting of the proton Ip levels in 
160. 

We have shown that the largest part of the spin-orbit effects in HF calculations is 
produced by the traditional shell model term, lAc in eq. (|^). The contribution of the 
other term, UlS) is very small and it can be simulated by a redefinition of the parameters 
of the force. The role of the exchange term is more complicated: its inclusion in the 
calculations modifies by a maximum of 15% the values of the spin-orbit splittings. The 
complication arises because these modifications do not have the same sign for all the nuclei 
studied. In calculations done with only the direct terms, it is not possible to simulate the 
exchange effects by simply readjusting the parameters of the interaction. 

Forcing our formalism to handle zero-range spin-orbit interactions we have studied, 
by comparison, the importance of the finite range. We found that calculations done with 
zero-range interaction produce energy splittings which, in heavy nuclei, are much larger 
than the empirical ones. Traditional HF calculations use spin-orbit zero-range terms 
of Skyrme type 0. These expressions produce contributions to the hamiltonian which 
are related to the derivative of the density distribution, while our expressions depend 
directly from the density distribution. The dependence from the derivative of the density 
distribution simulate effects produced by the finite range of the force. 

We have done calculations with spin-orbit terms taken from microscopic interactions 
[|, |§] and we have obtained splittings close to those produced by our effective spin- 
orbit interactions. This would indicate that the medium does not affect the spin-orbit 
term of the realistic interaction, in agreement with the findings of G-matrix calculations 
[0. We would like to point out, however, that the observables we have investigated 
are more sensitive to the global properties of the spin-orbit potential than to its details. 
Modifications of the local properties of the interaction would not produce effects on our 
results. 

We have also investigated the effects of the spin-orbit isospin term of the interaction, 
and we found them to be very small. These terms are related to the differences between 
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protons and neutrons density and spin-orbit density distributions. In our calculations 
these quantities are rather small even for a nucleus with a large neutron excess like ^''^Pb. 
There are however indications for observables which seems to be sensitive to this part of 
the potential . 



The comparison of the results of our calculations with the empirical values of the 
splittings on the various nuclei investigated is not satisfactory. The empirical splittings 
are smaller than those we have obtained, except for ^^C. This is a well known problem of 
the HF theory, and it could be solved only by using theories going beyond the mean field 
description of the nucleus. 



5 Appendix 

Since we have developed the HF equations in spherical coordinates it is necessary to 
express the operator L ■ S in terms of these coordinates. For this purpose we define an 
operator 0{ijk) as: 



^i-iy^^r, [V, ® aik)t^ , (21) 



where we have set h = 1 and the indexes i,j, k can assume only two values, 1 and 2 for 
example. The previous formula has been obtained by expressing r in terms of spherical 
harmonics and by making explicit use of the tensor product properties. Using the above 
operator we can write the spin-orbit channels of the force as: 

V;(ri2)Op(l,2) = Vp{n2) I" E 0{tjk) 

i,j,k=l,2 

^ ^ L jk 

J2{-lf'+'r,YL-M{ri)YLM{r2)Y,^{r,), p = 7,8. (22) 

iM 

The operator has been defined as: 

1, p = 7 

(23) 

r(l)-T(2), p = 8. 



In eq. (p^) we make the sum on M e /i, and we obtain a more synthetic expression: 



V;,(n,)0,(l,2) = -2v^X^i:f (n n' i )V£(ri,r2) 







EOo'o^(^jfc), P = 7,8, (24) 

ijk 
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where we have defined new set of operators as: 
and where we have defined: 

= 



ri, for i=2 
for i=l . 



(25) 



(26) 



It is convenient to express these operators such as the coordinates of each particle are 
separated: 

\L+M ILL' 1 [ 



Oor(uA:) = v^E(-l)' 



KM 



1 1 K 



(27) 



where we have used the Racah 6-j symbol and we have defined the operators: 



0^W(u^)=e(.jfc)ASl?(l)i?S(2). 



f28) 



The expressions of the three terms of the above equation, for each value of {ijk) are given 
in tab. as functions of the following operators: 



[YL{r^ 



(29) 



At the end we express the spin-orbit terms of the interaction as: 



LL'K 



ijk M 



where / is given by: 



(30) 



In the calculation of the HF equations for the spin-orbit channels we have used the 
results corresponding to the matrix elements 0^i^,j^{ijk) which in the tab. |^ are shown to 
be function of Cf^^, A^f^; ^lm defined in (p9|) and of the spherical harmonics Y^Mif)- 

Using the function ^(/) =1 if / is even and =0 if / is odd, we express the reduced matrix 
elements for the spherical harmonics as: 



{i\3\\yL\\i'\f) 



L 



-i 



(31) 
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For the other three operators we obtain the following operators: 



2^ 



'^{-rf"3 'j'lVLL'K { 



I' \ J' 
K 1 L 



+ 



/' K I 




K 1 L' 




K 1 L' 




dr I 



(i\j\\c^\\i'\j') 



47r 



■1) 



+ 



r L I 




L 1 X 




d ■ 
dr 



L 
1 



1 K 
-1 



^{-lyjfU'LK < 



I' I f 

1 L 



K 



I K V 




(32) 
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Table 1: Parameters of the central force 51a compared with the original parameterization 
of the Brink and Boeker Bl force [jll|. The values of the A are expressed in MeV, while 
those of b in fm~^. 





Bl 


Bla 


An 


-55.12 


-61.51 


Al2 


647.06 


647.06 


Api 


17.10 


15.82 


Ap2 


51.51 


51.51 


bl 


0.51 


0.51 


b2 


2.04 


2.04 



Table 2: Binding energies, in MeV, for the nuclei investigated. The force Bla does not 
have spin-orbit terms, the z, c, d and so calculations are discussed in the text. The 
interactions labeled U, A14 and A18 have been obtained by adding to Bla the spin-orbit 
terms of the Urbana V14 potential , of the Argonne V14 [^] potential and the Argonne 
VIS Q potential. The experimental energies are taken from ref. ||12|| . 







160 


40Ca 


^^Ca 


208pb 


Bla 


-48.89 


-116.09 


-341.54 


-380.48 


-1854.00 


z 


-54.47 


-116.08 


-341.43 


-402.67 


-2053.15 


c 


-55.26 


-116.09 


-341.52 


-397.74 


-1913.47 


d 


-63.13 


-116.42 


-342.30 


-413.36 


-1985.81 


so 


-67.37 


-116.42 


-342.24 


-420.60 


-2007.25 


U 


-62.77 


-116.46 


-342.27 


-411.75 


-1969.91 


AU 


-55.18 


-116.19 


-341.72 


-395.20 


-1907.95 


A18 


-61.20 


-116.37 


-342.10 


-408.55 


-1960.31 


exp 


-92.16 


-127.62 


-342.05 


-416.00 


-1636.45 
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Table 3: Spin-orbit splittings in MeV for the nuclei investigated. The empirical values have 
been calculated using the compilation of ref. [|14|. vr stands for protons, v for neutrons. 





z 


c 


d 


so 


U 


A14 


A18 


Sill 


exp 






















Ipvr 


3.37 


3.90 


4.97 


5.76 


4.88 


1.58 


4.16 




14.00 




3.05 


3.57 


4.65 


5.46 


4.61 


1.24 


3.83 




13.76 


160 




















Ipvr 


6.30 


6.31 


6.31 


6.32 


5.78 


2.86 


5.33 


5.90 


6.33 


Ipz/ 


6.38 


6.36 


6.36 


6.37 


5.89 


2.91 


5.41 


6.03 


6.16 



40Ca 




















Ipvr 


9.38 


5.56 


5.64 


5.35 


4.91 


2.40 


4.59 


3.31 




Idvr 


12.14 


8.36 


8.45 


9.16 


7.52 


3.70 


6.99 


6.18 


7.17 


Ipz/ 


9.47 


5.58 


5.66 


5.38 


5.03 


2.47 


4.67 


3.38 


3.00 


Idz/ 


12.43 


8.46 


8.56 


9.30 


7.74 


3.82 


7.15 


6.33 


6.30 


^«Ca 




















Ipvr 


10.05 


4.62 


5.11 


4.94 


4.47 


2.10 


4.13 


2.72 




Idvr 


13.97 


7.50 


8.27 


9.03 


7.52 


3.70 


6.99 


6.18 


4.30 


Ipz/ 


9.25 


3.89 


4.35 


3.76 


3.41 


1.23 


3.20 


2.36 




Idz/ 


12.68 


6.35 


7.09 


7.49 


5.98 


2.21 


5.47 


5.50 


3.60 


Ifz/ 


11.69 


7.82 


8.53 


10.17 


8.07 


3.28 


7.16 




8.74 
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Table 4: The same as tab. |for ^ospb. 





-y 




A 

LL 




TT 

yj 


A14 


AIR 


SIII 


CA.p 


2U8pu 






















IpTT 




i.yo 


9 dp; 


1 fin 

i .ou 


i .O i 


n P7 

U.O ( 


1 

i .o^ 


n 
U. 


c;9 
,0z 




1 Att 
iQTT 


99 ^^9 


0.00 


"X P'X 
0.00 


"X KP 
0.00 


9 7Q 

z. ( y 


i.Uo 


9 PP 
Z.OO 


i. 


91^ 
,Z0 




1 f-TT 

LL71 


98 ZL7 


O.Ui 




o .oo 


^.OO 


1 fi7 

i.O ( 


zL 1 n 


Z. 


Q/1 
.o4 




Igvr 


^^9 IP, 


p. Q7 

o.y ( 


7 77 


O.OO 


P OP 

o.zo 


z.oo 


O.OO 


Q 
O. 


, / 




1 Tttt 
illTT 




Q 1 
y. io 


1 n 1 7 

iU. i ( 


1 1 /17 












c;7 

0.0 ( 


ZpTT 


iZ.OO 






i.yo 


9 1 9 
Z. IZ 


n Qi 
u.y i 


1 QP 
i.yo 


n 


99 




2d7r 


18.78 


4.00 


4.18 


4.29 


3.60 


1.63 


3.36 


1, 


.76 


1.32 


2f7r 


21.03 


5.23 


5.46 


5.98 


4.77 


2.00 


4.41 


2. 


,35 


1.92 


Ipz/ 


14.63 


2.05 


2.11 


1.46 


1.80 


0.83 


1.68 


0. 


,53 




Idz/ 


22.44 


3.48 


3.74 


3.55 


3.17 


1.38 


2.94 


1. 


,26 




Ifz/ 


28.63 


5.10 


5.64 


5.83 


4.77 


2.00 


4.41 


2. 


,35 




Igz/ 


32.85 


6.96 


7.77 


8.38 


6.67 


2.79 


6.13 


3, 


.81 




Ihz/ 


34.06 


8.89 


9.92 


11.08 


8.79 


3.75 


8.00 


5, 


.56 




2pz/ 


12.50 


2.31 


2.43 


1.58 


2.08 


0.89 


1.93 


1. 


,08 




2dz/ 


18.74 


3.92 


4.10 


3.96 


3.73 


1.68 


3.40 


1, 


,98 




2fz/ 


20.94 


5.15 


5.36 


5.97 


5.17 


2.41 


4.64 


2. 


,78 


1.77 


3pz/ 


8.05 


2.01 


2.01 


1.52 


1.95 


0.92 


1.76 


1. 


,03 


0.94 



Table 5: Differences in MeV between the splittings presented in tabs. ^ and | and those 
obtained leaving out the isospin part of the spin-orbit interaction for Urbana, Argonne 
V14 and Argonne V18 potentials. We give here the minimum, maximum and average 
difference, in absolute value. 





Minimum 


Maximum 


Average 


160 


0.010 


0.030 


0.018 


40Ca 


0.000 


0.070 


0.032 


^^Ca 


0.020 


0.290 


0.088 


208pb 


0.000 


0.450 


0.128 
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Table 6: Expressions of the OLL'xi'^j^) operators. 



ijk 


eiijk) 






111 




CL'K 
^LM 


Yl-m 


112 


ri 


nL' 




121 


ri 




nL 


122 


^ V 


Yl'm 


cLK 
^L'-M 


211 


^ v 


CLK 
^L'M 


Yl'-M 


212 


r2 


nL 




221 






'^K-M 


222 


-•■4 




qL'K 
^L-M 
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Figure 1: Radial dependence of the spin-orbit interaction used. In the panel I our gaussian 
interactions are shown while in the other two panels the radial dependence of the spin-orbit 
terms (II) and spin-orbit isospin terms (III) of the Urbana |^ Argonne V14 |^ and Argonne 
V18 Q interactions. 



18 



